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Abstract.  Resonant  tunneling  is  employed  to  probe  the  inhomogeneous  strain  in  silicon-based 
quantum  dots.  When  submicron  structures  are  etched  from  a p-Si/SiGe/Si  double-barrier  het- 
erostructure, the  resonant  / ( V)  peaks  shift  and  develop  a fine  structure  consistent  with  pronounced 
strain  relaxation  in  the  SiGe  quantum  well.  We  calculate  the  strain  dependence  on  dot  size  by  finite 
element  techniques  and  convert  the  strain  to  an  effective  lateral  confining  potential.  In  sufficiently 
small  dots,  we  find  that  the  inhomogeneous  strain  confines  carriers  not  only  to  the  central  core,  as 
in  GaAs-based  dots,  but  also  to  a ring-like  region  at  the  perimeter.  We  probe  the  resulting  density 
of  states  by  magnetotunneling  1(V,  B ) measurements. 


Introduction 

When  a semiconductor  nanostructure,  like  a quantum  wire  or  a quantum  dot,  is  fabricated 
from  strained  epitaxially  grown  material,  the  originally  homogeneous  strain  is  replaced 
by  geometry-specific  strain  gradients.  The  symmetry-based  analytic  treatment  [i]  for 
handling  the  biaxial  strain  in  lattice-mismatched  materials  no  longer  suffices,  particularly 
for  quantum  dots  which  lack  translational  symmetry  altogether.  Instead  one  turns  to  finite- 
element  calculations  of  the  strain  field  based  on  linear  elastic  models,  but  their  applicability 
to  structures  whose  size  D might  be  down  to  tens  of  lattice  constants  is  not  self-evident. 
Thus,  any  experimental  technique  that  is  sensitive  to  inhomogeneous  strains  provides  a 
valuable  test  bed  for  the  validity  of  finite-element  techniques  for  strained  nanostructures  — 
in  our  work,  the  experimental  probe  will  be  the  resonant  tunneling  current- voltage  I(V) 
measurements. 

Further,  inhomogeneous  strain  in  semiconductor  nanostructures  is  taking  on  additional 
technological  relevance,  as  advances  in  strained  layer  epitaxy  and  ongoing  device  miniatur- 
ization promise  the  arrival  of  deep  submicron  bandgap-engineered  devices,  such  as  strained 
Si/SiGe  HBTs  for  high-frequency  analog  and  digital  applications.  Particularly  interesting 
is  the  use  of  strain-driven  self-assembly  of  quantum  dots  in  semiconductor  lasers  [ J to 
enhance  gain  and  shift  the  lasing  wavelength.  In  such  devices,  strain  relaxation  is  the  key 
to  the  size  and  morphology  of  the  quantum  dots,  as  well  as  the  optical  transition  energies. 

Finally,  it  is  important  to  note  that  inhomogeneous  strains  in  quantum  nanostructures 
can  contribute  to  carrier  localization  in  unpredictable  ways.  Over  the  past  decade,  quantum 
dots  have  been  extensively  investigated  as  systems  containing  a few  spatially  confined 
charge  carriers  [3-5].  However,  most  of  these  experiments  probed  dots  made  from  lattice- 
matched  GaAs/AlGaAs  heterostructures,  in  which  the  carriers  are  confined  to  the  central 
region  of  the  dot  by  a roughly  parabolic  lateral  potential  arising  from  the  gate  potential  or 
the  pinning  of  the  Fermi  level  at  the  surface.  In  our  strained  Si/SiGe  quantum  dots,  the 
inhomogeneous-strain-induced  lateral  confinement  potential  is  nonmonotonic,  leading  to 
an  effective  potential  minimum  near  the  perimeter  of  the  dot.  For  sufficiently  small  dots. 
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Bias  (mV) 

Fig.  1.  (a)  SEM  micrograph  of  a representative  double-barrier  nanostructure  with  lateral  dimension 
D ~ 0. 15  pm.  (b)  Self-consistent  potential  profile  of  the  Si/SiGe  double-barrier  active  region  at 
V = 150  mV.  (c)  Resonant  tunneling  /(V)  of  a large  D = 2 /urn  device  at  1 .7  K,  with  the  HH  peak 
magnified  x 35  for  clarity. 

this  potential  confines  the  ground  state  to  a one-dimensional  (ID)  ring-like  region  near  the 
perimeter  — a new  type  of  structure  expected  to  have  interesting  magnetic  properties. 

1.  Tunneling  and  strain  in  SiGe  double-barrier  microstructures 

Our  devices  begin  with  p-Si/Si]_xGeA/Si  double-barrier  heterostructures  described  in  detail 
in  previous  publications  They  are  grown  on  p- Si  substrates,  with  an  undoped 

active  region  consisting  of  Si  barriers  confining  a Si|_vGex  QW  that  is  35  A wide  with  Ge 
content  x = 0.25  or  0.2  (corresponding  to  a lattice  mismatch  of  ~ 1 and  ~0.8%  respectively). 
Outside  the  barriers  are  p-Si|_AGe.Y  emitter  and  collector  regions  that  serve  as  reservoirs 
for  tunneling  holes. 

When  a bias  V is  applied  between  the  emitter  and  collector,  the  holes  in  the  emitter  tunnel 
via  the  quantized  2D  hole  subbands,  subject  to  the  usual  energy  E and  transverse  momentum 
kj  conservation  rules  [ ] . In  large  devices,  the  strain  in  the  SiGe  well  can  be  taken  as  biaxial 
and  homogeneous,  so  the  energies  of  the  2D  subbands  can  be  reliably  calculated  numerically 
| 1 . Figure  1 shows  an  SEM  photograph  of  a device,  together  with  a self-consistent  potential 
distribution  in  the  active  region  under  bias,  and  the  I (V)  characteristic  of  a large  D — 2 pm 
Si/Sio.75Geo.25  device  at  T — 1.7K.  The  peaks  correspond  to  tunneling  through  the  two 
confined  2D  subbands,  labeled  HH  and  LH  for  the  heavy-hole  and  light-hole  branches  of 
the  dispersion.  The  agreement  with  the  predicted  peak  positions  is  excellent  | , ],  so  the 
peak  positions  reflect  the  energies  of  the  quantized  states  in  the  SiGe  well.  In  particular, 
the  energy  separation  A E between  the  HH  and  LH  subbands  arises  in  part  from  the  strain- 
induced  splitting,  which  lifts  the  HH-LH  degeneracy  in  the  SiGe  valence  band  | , ]. 

Thus,  any  significant  change  in  the  strain  should  be  reflected  in  A E. 
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Fig.  2.  I (V)  characteristics  with  nominal  lateral  diameters  of  D = 1 .0—0.25  [im  at  1 .7  K.  Current 
scale  corresponds  to  the  smallest  device,  other  curves  are  rescaled  and  the  HH  peaks  magnified  for 
clarity.  Vertical  lines  mark  the  HH  and  LH  peak  positions  in  large  devices. 

2.  Average  strain  relaxation  in  SiGe  double-barrier  niierostructures 

When  a narrow  mesa  is  etched  through  the  active  region  of  our  devices,  the  biaxially 
compressed  SiGe  strained  layers  can  relax  by  lateral  displacement  at  the  sidewall.  We  first 
observed  this  effect  experimentally  some  years  ago  | ] and  it  has  been  confirmed  by  others 
| ].  Figure  2 shows  the  effect  in  a series  of  identically  fabricated  Si/Sio.75Geo.25  devices 

with  D ranging  from  1 down  to  0.25  ^un.  The  I(V)  curves  exhibit  a consistent  shift  of 
HH  and  LH  peaks  towards  each  other  as  D decreases,  even  as  the  I (V)  lineshape  remains 
largely  unaffected.  The  effect  is  large,  with  changes  in  A E indicating  a significant  amount 
of  strain  relaxation. 

We  compared  the  strain  relaxation  inferred  from  the  data  in  Fig.  2 to  finite-element 
simulations  based  on  a linear  elastic  model,  in  which  the  cylindrical  structure  was  allowed 
to  relax  to  a minimum  energy  configuration  [ ].  The  average  values  of  strain  relaxation 

predicted  by  these  calculations  agree  quite  well  with  our  experimental  data,  with  the  domi- 
nant radial  strain  component  €rr  relaxing  to  ~ 0.7  of  the  full  lattice-mismatch  strain  when 
D falls  to  0.3  n m.  Interestingly,  the  simulated  strain  relaxation  in  the  SiGe  layers  is  non- 
monotonic in  the  radial  direction  rj_  and  significant  strain  gradients  exist  throughout  the 
SiGe  well  in  sufficiently  small  structures,  D < 0.25  //m.  The  inhomogeneous  strain  leads 
to  additional  lateral  quantization  in  the  confined  2D  subbands,  which  we  first  observed 
in  | |. 

3.  Inhomogeneous  strain  in  SiGe  quantum  dots 

The  calculated  radial  strain  err  in  the  SiGe  QW  and  the  corresponding  strain-induced  lateral 
potential  for  HH  states  are  shown  in  Fig.  3 as  a function  of  r±  for  D — 0. 1 —0.2  //  m devices 
with  a Sio.8Geo.2  well.  The  inset  shows  the  calculation  geometry  and  lateral  displacement 
of  the  strained  SiGe  layers.  First  consider  the  radial  strain  err(r±)  curves  at  the  top  of 
Fig.  3.  For  D — 0.2  /<m,  err  decreases  gradually  with  r±  except  for  a region  of  increasing 
strain  near  the  perimeter  that  extends  ~ 100  A and  reaches  err  ~ 0.6. 

This  strained  ring-like  region  exists  for  all  D [ ],  [ ],  but  for  devices  much  larger  than 
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In-plane  radius  (A) 

Fig.  3.  The  top  curves  show  the  calculated  radial  strain  component  err(r±)  for  D = 0.1—  0.2  n m 
devices  at  the  mid-plane  of  the  Sio.sGeo.2  well  (full  strain  is  err  = 1).  Inset  shows  the  magnified 
displacement  of  the  finite-element  mesh  near  the  sidewall.  The  bottom  curves  are  the  corresponding 
in-plane  confining  potentials  for  the  HH  states  as  function  of  r±.  Dashed  lines  mark  the  confined 
ring  subbands  in  a D ~ 0. 15  jiim  dot. 

D — 0.2  /I m it  can  be  taken  as  a perturbation  to  the  inner  core.  On  the  other  hand,  in  small 
devices,  D < 0. 15/ttm,  the  ring  at  the  perimeter  becomes  the  most  highly  strained  region 
in  the  device.  The  corresponding  confining  potential  for  HH  states  is  shown  by  the  bottom 
set  of  curves  in  Fig.  3:  in  the  smallest  dots,  the  strain  would  confine  holes  to  a ID  ring  at 
the  perimeter. 

Figure  4 shows  the  HH  I(V)  peak  of  a D ~ 0.1 5 ,um  device  fabricated  from  the 
Si/Sio.8Geo.2  double-barrier  material,  together  with  a reference  lineshape  from  a large 
device.  The  I(V)  lineshape  exhibits  very  strong  fine  structure,  corresponding  to  strong 
lateral  quantization  due  to  inhomogeneous  strain  | ].  A full-blown  calculation  of  the 
expected  density  of  states  in  such  a quantum  dot  in  the  presence  of  the  inhomogeneous 
strain  is  complicated  by  the  anisotropy  and  nonparabolicity  of  the  in-plane  effective  mass 
in  the  quantum  well,  but  taking  the  in-plane  HH  effective  mass  m*  ~ 0.25  we  obtain  several 
radially  quantized  subbands  in  the  perimeter  ring  separated  by  a few  meV,  see  Fig.  3.  Since 
quantized  ring-like  subbands  overlap  in  energy  with  the  states  in  the  relaxed  central  core, 
they  would  be  expected  to  contribute  structure  on  top  of  the  relatively  smooth  overall  HH 
peak  lineshape.  The  energy  separation  of  the  ring  subbands  extracted  from  the  structure  in 
the  I(V)  agrees  reasonably  with  the  calculations  of  Fig.  3. 

4.  Magnetotunneling  spectroscopy  of  strained  SiGe  quantum  dots 

The  additional  ring-like  confinement  of  hole  states  in  sufficiently  small  inhomogeneously 
strained  SiGe  quantum  dots  provides  a new  and  interesting  system  for  magnetic  field  effects. 
Figure  5(a)  shows  the  evolution  of  the  HH  I (V)  peak  fine  structure  in  magnetic  fields 
B — 0— 10  T (B  ||  /),  while  Fig.  5(b)  indicates  the  evolution  of  the  I(V)  peaks  converted 
to  an  energy  scale.  In  the  absence  of  inhomogeneous  strain,  the  B field  would  compress 
the  2D  subband  density  of  states  into  Landau  levels.  Given  a constant  in-plane  effective 
mass,  as  in  u-GaAs/AlGaAs  double-barrier  structures,  the  Landau  level  separation  would 
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Fig.  4.  HH  resonant  peak  1(V ) and  d//dV  characteristics  of  a D ~ 0.15/xm  Si/Sio.sGeo.2  dot. 
The  dashed  line  shows  the  smooth  HH  I (V)  lineshape  in  a large  device. 


Bias  (mV)  B (T) 


Fig.  5.  (a)  Evolution  of  HH  I {V,  B)  lineshape  of  the  D ~ 0. 1 5 fx  m dot  with  B ; (b)  5-field-induced 
shifts  and  splittings  in  the  density  of  states  peaks  in  the  dot. 


be  linear  in  B and  the  / (V/,  B ) lineshape  would  show  equally  spaced  features  falling  on  the 
usual  Landau  fan  diagram  [ ].  Given  the  complex,  nonparabolic  dispersion  of  holes  in 

SiGe  quantum  wells,  the  Landau  level  spectrum  is  quite  complicated  even  in  a uniformly 
strained  well  [ | and  /(-induced  structure  has  only  been  seen  experimentally  at  fairly  large 

B | ].  The  data  in  Fig.  5 shows  rather  complex  behavior,  with  some  of  the  peaks  shifting 
towards  lower  energy  with  B,  others  appearing  to  repel  each  other,  and  some  even  splitting 
around  B ~ 5 T.  A theoretical  analysis  of  ring-like  hole  states  in  a B field  remains  to  be 
performed. 


5.  Conclusions 

We  have  investigated  the  effects  of  size-induced  strain  relaxation  in  strained  SiGe  quantum 
dots.  Our  data  indicate  that  large  strain  relaxation  and  nonomonotonic  strain  gradients  ap- 
pear in  deep  submicron  structures,  in  agreement  with  finite  element  simulations.  We  also 
see  evidence  for  confinement  of  carriers  to  ring-like  regions  at  the  perimeter  due  to  inhomo- 
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geneous  strain.  Our  measurements  prove  resonant  tunneling  to  be  a viable  spectroscopic 
probe  for  strain  effects  in  individual  nanostructures. 
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